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Abstract 

In the paper a possibility of metallic ferromagnetic state realization in a gen- 
eralized Hubbard model with more complete accounting of electron-electron inter- 
actions, in particular, the correlated hopping and exchange interaction integrals 
is investigated. Recently obtained by means of mean-field approximation single 
electron energy spectrum is used for the description of finite temperature prop- 
erties of the system. In the paper the expression for the critical temperature of 
ferromagnet-paramagnet transition is found, the behaviour of temperature depen- 
dencies of magnetization and paramagnetic susceptibility is analyzed. Taking into 
account the correlated hopping allows to explain some peculiarity of ferromagnetic 
behaviour of transition metals, their alloys and compounds. 

PACS number(s): 71.10.Fd, 71.30.+h, 71.27.+a 

1 Introduction 

The important puzzle in the explanation of ferromagnetism in the system of single band 
electrons is proper taking into account the Coulomb correlation between electrons. The 
Hubbard model [I], 0, which describes the single non-degenerate electron band with 
local Coulomb interaction, is incapable to describe ferromagnetism and requires the gener- 
alization. It is naturally to generalize the model Hamiltonian by accounting other matrix 
elements of electron-electron interaction (in addition to the intra-atomic Coulomb repul- 
sion) and to consider the ferromagnetism in the generalized Hubbard model (for review 
on this problem see Refs. ||, ||, |(|). 

Note that last time the problem of metallic ferromagnetism in the single band Hubbard 
model and their generalizations has attracted much attention of researchers in a series 
of papers by means of the Gutzwiller variational wave function approximation [0, ||, 
the dynamical mean-field theory ||, 10], the spectral density approximation [11|, the ex- 



act diagonalization method fl3| and mean-field theory |14 , |15|, [TJ|. In the papers of 



Hirsch |fH] , pT5| , |I"6|j the generalization of Stoner-Wolfarth theory [17], pl| was carried out 



The Stoner-Wolfarth theory has been used for the description of band electron magnetism 
(denote in this connection, that the calculations in this theory essentially depend on the 
shape and peculiarities of the density of states, in particular, it is known that the incom- 
plete ferromagnetism in this model is absent if the density of state is rectangular |18[). The 
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considered generalization [14, 15, 16] has shown that the inter- atomic exchange interac- 
tion plays the key role for the obtaining of ferromagnetic state with partially polarization. 
Using the mean-field theory Hirsch has obtained the condition of ferromagnetic state real- 
ization, nonmonotonic behaviour of the concentration dependence of magnetization, the 
expression for the Curie temperature, temperature dependencies of magnetization and 
magnetic susceptibility. The authors of Ref. || using local approximation developed from 
the Gutzviller wave function method also suggest (in agreement with |14]]) the importance 
of inter-atomic exchange interaction for the stabilization of ferromagnetic state with par- 



tial spin polarization. In the Ref. [19] the opposite case of strong correlation has been 



studied and the criterion of ferromagnetism and magnetization of the system in ground 
state have been derived. 

At our point of view it is in principle necessary to consider in addition to the inter- 
atomic exchange interaction the matrix elements that describe the correlated hopping of 
electrons |25| for the explanation of ferromagnetism in single band model. Taking into 
account the correlated hoppind has alloved to explain some peculiarities of metal-insulator 



transition in single band model [EG, 21] and double orbitally degenerate model [22]. In this 



model an electron hopping from one site to another is correlated both by the occupation of 
the sites involved in the hopping process and the occupation of the nearest-neighbour sites. 



Taking into account the correlated hopping allows to obtain the additional mechanism p3 



2|| for the stabilization of ferromagnetic ordering. 

The importance of correlated hopping for understanding of the metallic ferromag- 
netism was discussed in Ref. |Tj|, where the results, obtained by means of exact diagonal- 
ization and mean-field theory for small one-dimensional chain of atoms were compared. 
It has been shown that in strong coupling regime correlated hopping favours ferromag- 
netism stronger for electron concentration n > 1 than for n < 1 (the reverse situation 



occurs at weak interactions [131); this result agrees with conclusions of Ref. p5| . Using 
the Gutzviller approach the authors of Ref. J/J have shown, that the correlated hopping 
favour the ferromagnetic ordering around the point of half-filling. 

The magnetic properties of the system at nonzero temperature in the case of strong 



and intermediate electron correlations was analyzed in a series of papers [§, [L0| [11], |12]| . 
However the important question concerning the influence of correlated hopping on the 
Curie temperature, on the behaviour of magnetization and magnetic susceptibility at 
nonzero temperature is not yet considered. 

The main goal of the article is to show the principle necessity to take into account 
both correlated hopping of electrons and inter-atomic exchange for the correct description 
of metallic ferromagnetism. Note that the application of mean-field theory leads to the 
overestimation of ferromagnetic tendencies and gives the results, which agree with exper- 
iment only qualitatively. At the same time it is interesting to investigate the case of weak 
electron correlations and to define the mechanismes which are responsible for the ferro- 
magnetic ordering. The mean-field analyzis in spite of its limitation allows to reproduce 
the behaviour of the magnetic moment and Curie temperature at the change of electron 
concentration, the correct temperature dependence of magnetic moment of the system. 
Besides we believe that taking into account the correlated hopping allows to avoid the 
overestimation of Curie temperature. 

The paper has the following structure. In the Sec. 2 the Hamiltonian of the generalized 
Hubbard model with correlated hopping and inter-atomic exchange interaction is written 
and the single electron energy spectrum obtained in the case weak electron interactions by 
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means of Green function technique is analyzed. In the Sec. 3 the ground state properties of 
the system is considered, the equation on the system critical parameters and expression for 
the magnetization are derived. In the Sec. 4 the finite temperature properties of the system 
are considered, the expression for the Curie temperature, temperature dependencies of 
magnetization and magnetic susceptibility are obtained. The Sec. 5 is devoted to the 
conclusions. 



2 The model Hamiltonian and single electron en- 
ergy spectrum 



Let us write the Hamiltonian, proposed in Refs. |p0| , |23|j , generalized by taking into account 
weak magnetic field 

H = - n a ta a i° + 12 t ij( n ) a ta a j^ + 12' {T 2 {ij)at a a ja n ia + /i.e.) 

ia ija ija 

+ rkfnu + ~Y1 'af a a+ a ,a ia ,a j(T - (n^ - n^) , (1) 

where af a , a ia are creation and destruction operators of electron on site i, a =|, I, n ia = 
af a a ia , n ia = af a a ia . is the operator of number of electrons with spin a on site i, n = 
(n^+riii), /i is the chemical potential, £y(n) = tij+nTi(ij) is the effective hopping integral 
of an electron from site j to site i, tij is the band hopping integral an electron from site j 
to site i, T\{ij) and T 2 (ij) is the parameters of correlated hopping of electrons, U is the 
intra-atomic Coulomb repulsion, J is the exchange integral for the nearest neighbours and 
h is external magnetic field (the units of h use such the magnetic moment per electron 
is unity). The prime on the sums in Eq. (JJ) signifies that i ^ j. The concentration 
dependence of the effective hopping integral (n) is caused by the correlated hopping p(J 
of electrons. 

The peculiarities of the model described by the Hamiltonian (|l|) are taking into account 
the influence of site occupation on hopping process (correlated hopping) and the direct 
exchange interaction between electrons in neighbour sites. In this model an electron 
hopping from one site to another is correlated both by the occupation of the sites involved 
in the hopping process and the occupation of the nearest-neighbour sites (in contrast to 
similar generalized Hubbard models). 

To characterize the value of the correlated hopping we introduce dimensionless param- 



Tijij) ^ = T 2 (ij) 

In terms of Green functions method using the mean-field approximation we have re- 



efers T\ = }}y , T2 = f, i which are independent on the number of site. 



cently obtained |2J] the single electron energy spectrum, which in the external magnetic 



field has the following form 

E a (k) = —[i + (3 a + n„U - zn a J - hr] a + t(na)7(k); (2) 

here the spin-dependent shift of the band centers is 
2 

& = ]y 12 T 2 (y ) ( a i* a 3°) > ( 3 ) 

ij 
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z is the number of nearest neighbours to a site, for the spin o =| (!) we have t] a = 1 (— 1), 

7(k) = £e ikR (the sum goes over the nearest neighbours to a site), and the spin and 
R 

concentration dependent hopping integral is 

zJ 

t(na) = (1 - T%n - 2r 2 n s ^ 4r0* = ( 4 ) 

^ cr' 

where w = z\t\ is the half of band width, t is the hopping integrals between nearest 
neighbour sites and 

A, = ^E(-^<^v>). (5) 

The energy spectrum (§) will be used in the next sections for the description of the model 
properties in the ground state and at the finite temperature. 

3 The ground state properties of the model 

In general the occupation number and the magnetization are expressed respectively as 

w 

n = n T + n { = I p(e) [/(_&[■ (e)) + /(^(e))] de; (6) 



m = n T - ni = J p(e) [/(£ T (e)) - /(^(e))] de, (7) 

—w 

where p(e) is the density of states, f(E a (e)) is the Fermi function. 

Generalizing the recently obtained result |24| (for the rectangular density of states) in 
the case of weak external magnetic field one can obtain for the correlation function A a 
and the shift of the center of u-subband f3 a 

A a = n a (l-n a ), (8) 

(3 a = 2wr 2 A a = 2wT 2 n a (l - n a ), (9) 

where the concentration of electrons with spin a 

£* + w 

n a = (10) 

here e a = ^- is the solution of the equation E a (e) = 0, where p a = p — j3 a + zn a J — n & U + 
hr] a and a a = 1 - r x n - 2r 2 n s - ^ E A/- 

a' 

On basis of expressions (§|)-(0) using (§)-(§) one can obtain the equation on the system 
critical parameters 

Z A (i + n{2 - n )-m 2 )=l- nr x - r 2 (2 - n) - ^- - -. (11) 
2w v ' 2w w 

The condition for onset of ferromagnetic ordering is obtained (replacing the sign of equality 
by sign of inequality), when m — 0, and the condition for the ferromagnetic state with 
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full spin polarization if one put m = n in the expression (|TT|). The peculiar distinction 
of the obtained expression from the similar condition, obtained in Refs. |TJ], [L5], [TJ^ |8| for 
the various generalizations of the Hubbard model is the presence of correlated hopping, 
which can essentially modify the properties of the model. Taking into account correlated 
hopping leads to the appearance of a peculiar kinetic mechanism of ferromagnetic ordering 
stabilization. This mechanism is caused by the presence of the spin-dependent shift of the 
spin subband centers being the consequence of correlated hopping (which are analogous to 
the shift of subband centers in consequence of inter-atomic direct exchange interaction). 

The spontaneous magnetization (in the absence of magnetic field) of the system is 
found from the equation flTlD as 



m GS 



„ (2 _ B) _ .-^-,»-.(2-n) [ (12) 



which is valid in the case J > (when J = only the transition from paramagnetic to 
fully polarized ferromagnetic state with m GS = n occurs). If the calculated magnetization 
m GS > n it is necesarry to put m GS = n. 

The influence of correlated hopping on the properties of the system is illustrated on the 
Figs. 1, 2. On the Fig. 1 the dependence of the critical value of exchange integral at which 
the ferromagnetic ordering occurs is plotted (at the absence of magnetic field) as a function 
of band filling at various values of correlated hopping parameters r 1; r 2 and intra-atomic 
Coulomb interaction U/w. Solid curves correspond to the onset of spin polarization, 
dashed curves correspond to the fully polarized ferromagnetic state (the area below the 
solid line is paramagnetic, above the dashed line - fully polarized ferromagnetic, between 
lines - partially polarized ferromagnetic). The similar phase diagrams have been obtained 
in the Refs. |14] , [16| , but in those works the correlated hopping was not considered. Note 
that the presence of correlated hopping (curves 2, 3 on the Fig. 1) leads to the essential 
change of phase boundary of paramagnet-ferromagnet transition, in particular, to the 
shift of the minimum point, namely to the inequivalence of the cases n < 1 and n > 1. 

On the Fig. 2 the dependence of ground state magnetization m GS on the electron 
concentration n at the fixed values of exchange and Coulomb interactions, as well as 



at different values of correlated hopping parameters is plotted. In Refs. |T4], [L6| the 
concentration dependencies of magnetization in the framework of Hubbard model with 
inter-atomic exchange were obtained, but they are symmetrical respecting the point of 
half-filling. In the Ref. fl9j (for the case U — > oo) the similar m GS {n) dependence was 
obtained. 

As mentioned above, the asymmetry of cases n < 1 and n > 1 is obsereved; in 
particular, the increase of parameter n leads to the shift of the ferromagnetic area in 
a region of larger electron concentration n, the increase of parameter r 2 - in a region of 
smaller n. Note also that taking into account the correlated hopping considerably enriches 
the set of curves m GS (n). The obtained concentration dependencies of magnetization 
allow qualitatively describe the experimental curves for the binary ferromagnetic alloys 
of transition metals Fe, Co, Ni (Slater-Pouling's curves [f27|j ). 
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4 The properties of the model at nonzero tempera- 
ture 



At nonzero temperature and rectangular density of states the concentration of electrons 
with fixed spin direction is expressed as 



1 



e 



-In 



1 + exp(E a (w)) 



2wa a \1 + exp(E cr (—w)) J 
Using this expression one can obtain the equation for the magnetization 

mJ e f f w\ _ sn { e ) 311 \~~e~) 



(13) 



exp 







e 



(14) 



where J e // = zJ/w + U/w + 2r 2 (l — n). To obtain the temperature dependence of mag- 
netization it is necessary to apply the numeric methods forasmuch the last equation can 
not be solved analitically. The numeric calculations show that the results could be ap- 
proximately expressed using the approach proposed in Ref . [|15| . The critical temperature 
of transition from paramagnetic to ferromagnetic metallic state (the Curie temperature) 
can be obtained by expanding (0) at lowest order of m — * 0. 



1 U + zJ - 2wT 2 (n - 1) 



2w 



a" 



dfjx - fi* 
dx 



dx 



r-2 



w(a 



*\2 



here x = a*e, a* 



dfjx - fi* 
dx 



= o and fi* = fieri 



rn=0 • 



(15) 



At low temperature one can approximately write 

U/2w + zJ/2w - 2wr 2 {n - 1) 
~ 1 - (n + r 2 )n - 2A*(Q c )zJ/w' 



(16) 



where A*(Qc) — A cr (Q)\ m= o, = &bT, ks is Boltzmann constant. In the case of zero 
temperature the last expression reproduces the criterion of the paramagnet-ferromagnet 
transition, obtained in the recent paper p4| . 

To find Aer(Q) one can apply the Sommerfeld expansion |26| (at low temperature) 



AeriQ) 



2w 2 



de 



exp(E a (e)/Q) 



(17) 



71 

~3 



1 ■ 0\ 2 
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In the Curie point (when the magnetic moment m 
is written as 



2wJ [1 - Tl n - 2m- a - (A T (0) + A^Q^zJ/w} 2 ' 

and n a = n/2) the last equation 



A*(e c ) 



n(2 



n 



2ir 2 /e, - 2 



3 \2wJ [l-(T 1 + T 2 )n-2A*{e c )zJ/wY 
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Solving the system of equation (|T6|) and ([18|) we can express the Curie temperature as 
function of model parameters 



6c l-UI2 W -nn-r 2 (2-n) 
2w \ zJ/2w 



i^ + £^-2r 2 (n-l)). 119) 



V 2tt 2 V2w 2w 

Taking into account the expression ([12] ) we finally obtain for the Curie temperature 



2m V2 ir n2u, + ^- 2T2(n - 1) J m (20) 

At zero correlated hopping the formula ( pOf ) reproduces the result of Ref. [15|]. Note 



that the Curie temperature is tightly related (in the used approximation) to ground state 
magnetization of ferromagnet. The peculiarity of formula fl20|) is the presence of term 
related to correlated hopping r 2 (m GS is also dependent on this parameter), which in the 
case n > 1 is negative and therefore lowers the Curie temperature. This circumstance (at 
correspondent choice of correlated hopping parameters), at our opinion allow to overcome 
the considerable defect of mean-field treatment, namely the overestimation of the Curie 
temperature (in this connection also see Ref. ]nj, where the correlated hopping is not 



took into account). Let us put (as in the Ref. [jn|) the bandwidth of the e g states in 
3d ferromagnetic transition metals is approximately 2 eV. In our work the band filling 
n = 1.2 corresponds to Fe and for the values of correlated hopping n = 0.15, r 2 = 0.2 
our theory would predict the Curie temperature between 1000 and 1600 K depending on 
the values of U/w and zJ/w (these parameters are changed in the region from to 0.4 
and from 0.37 to 0.14 respectively). The value of Curie temperature at some fixed values 
of intra-atomic Coulomb interaction and inter-atomic exchange interaction agrees with 
experimental data even quantitatively. 

On the Fig. 3 the concentration dependence of the Curie temperature is plotted at 
various values of correlated hopping. The peculiarity of this dependence is the lowering 
of the Curie temperature at the increase of carriers concentration. Besides taking into 
account the correlated hopping causes the asymmetry of curve respecting the point of 
half-filling, that allows to explain qualitatively the larger Curie temperature in Co as 
compared to Fe. Note for the comparison that this fact can not be explained in the 
Ref. jL5) without further comments (in particular, taking into account the density of 
states with peculiarities). On basis of obtained expression for the Curie temperature 
one can explain the peculiarities of Curie temperature behaviour in the binary alloys of 
transition metals p7| . It is reasonable in the framework of our theory to interpret the 



paramagnet-ferromagnet transition in metallic phase at the increase of temperature for the 
nonstoichiometric chalcogenide chrome spinel Cro.5Feo.5S (where the Curie temperature 
is of the order of 1000 K) g§ . 

In principle, our theory allows to obtain the concentration dependences of magne- 
tization and Curie temperature, which are similar to the experimentally observed in 
compounds Fei-^Coa^ and Coi^Ni^^ with the change of electron concentration in 
3<i-band [p5fl . m these crystals the same subsystem of electrons is responsible both for 
conductivity and for the localized magnetic moments formation. Although these com- 
pounds are adequately described in the framework of double orbitally degenerate model, 
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even in terms of single band model at some values of model parameters it is possible to 
obtain the concentration dependencies of the Curie temperture and magnetic moment 
of system (see. Fig. 4) that agree qualitatively with experimentally observed f29fl . The 



plotted curves show that at choosen model parameters the values of mentioned above 
quantities are reproduced accurately. 

The influence of model parameters on the critical temperature is illustrated in the 
next figures. On the Figs. 5 and 6 the dependencies of the Curie temperature on the 
parameter of intra-atomic Coulomb interaction are plotted at half-filling (Fig. 5) and 
various values of band filling (Fig. 6), the values of exchange integral and correlated 
hopping are fixed. The plotted curves have a peculiarities: one can distinquish the area 
of sharp increase of the Curie temperature at the increase of parameter U/w (this values 
of model parameters correspond to the partial spin polarization of system) and the area, 
where the Curie temperature changes in proportion to U/w (this values correspond to full 
spin polarization of system). Note that the increase of inter-atomic exchange interaction 
leads to the extension of partially polarized ferromagnetic area, and also to the decrease 
of critical value of U, required for the development of saturated ferromagnetism. The 
peculiarity of the dependence of the Curie temperature on the U/w is the increase of 
critical value U / w at the increase of n. Note also that the changing of this critical value 
(at choosen values of correlated hopping) in the case n > 1 and at the increase of electron 
concentration is more pronounced than in the case n < 1. It is also interesting that 
depending on the magnitude of U/w the system with electron concentration n < 1 can 
have the larger value of the Curie temperature than the system with n > 1. However 
starting from the critical point for the intra-atomic Coulomb interaction the situation 
becomes opposite: the Curie temperature is larger for the system with n > 1. 

Next consider the behaviour of the magnetization of the system at the change of 
temperature. To obtain the temperature dependence of magnetization m let us use the 
assumption (as in the Ref. HT5[), that the equation (^) is also valid for the nonzero 
magnetization and for the temperature lower than critical Be- Than using (|T7|) one can 
obtain the equaiton 



2tt 

m 



2 1 1 / Qn2 

(21) 



e c \ 2 i we\ 



\2w) (a*) 2 2 \2wJ 



1 1 



(a* + T2m) 2 (a* — T2m) 2 



The results of numeric calculations for m is plotted on the Fig. 7. It should be noted that 
if the correlated hopping is not take into account, the last expression gives the analitical 
result of Ref. [|l^]. In the band limit last expression reproduces the result of Ref. [00(1 , 



where the similar treatment is applied to the completely itinerant carriers. 

To find the magnetic susceptibility let us to take the deriviative from magnetization 
(|7|) with respect to the magnetic field 

dTn(li) 1 1 

X(T) = ~0/T l '~ o ' m -* = 2w~ a* - U/2w - zJ/2w + 2r (n - 1) ' (22) 



Using ( |16|) at temperatures close to Oc one can obtain 
_ 6w[U/2w + zJ/2w-2r 2 (n-l)] 2 

x{1) ~ 27r 2 e c (e -e c )zj/w ' 1 ' 

On the Fig. 7 the temperature dependencies of magnetization and inverse magnetic suscep- 
tibility are plotted. The similar plot was obtained for the transition metals in the Ref. [ETT| 
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using the dynamical mean field theory with local density approximation. At temperatures 
higher than Curie temperature the magnetic susceptibility demonstrates Curie- Weiss-like 
behaviour, intrinsic for the ferromagnetic 3d-metals and their alloys. The similar temper- 
ature dependence of magnetic moment in alloy MnSb x As\- x at x > 0.07 and the presence 



of magnetic field is observed [32 



5 Conclusions 

In this paper the ferromagnetic solution in a single band generalized Hubbard model with 
correlated hopping is derived. The peculiarity of the model is the consideration so-called 
off-diagonal matrix elements of Coulomb interaction between electrons, in particular the 
integrals of correlated hopping (which describe the influence of occupancy of sites on the 
hopping process) and inter-atomic exchange interaction. The physical mechanisms that 
lead to the ferromagnetic state realization are the band narrowing and shift of the spin 
subband centers, caused by the exchange interaction and correlated hopping. 

Taking into account the correlated hopping leads to the asymmetry of cases n < 
1 and n > 1 at the consideration the ferromagnetism in this model. The increase of 
parameter T\ leads to the shift of the ferromagnetic area in a region of larger electron 
concentration n, the increase of parameter r 2 - in a region of smaller electron concentration 
n. The important consequence of the performed study is the conclusion that for the 
ferromagnetism realisation the case of more than half-filled band is more favourable than 
less than half-filled one;this maps the situation transition metals and their alloys. 

The concentration dependence of magnetization m GS qualitatively agrees with exper- 
imental data for the transition 3d-metals and their alloys, in particular, our results can 



explain the Slater-Pouling's curves p7| for the binary ferromagnetic alloys of transition 
metals Fe, Co, Ni with other 3d-metals. 

Taking into account the correlated hopping at the calculation of the Curie temperature 
allows to avoid the problem of the overestimation of critical temperature which is typical 
for the standard mean-field treatmens. The calculated Curie temperature is character- 
ized by peculiarities of concentration dependence, in particular, the asymmetry (in conse- 
quance of correlated hopping accounting) respecting to the point of half-filling. This fact 
allows to explain the experimentally observed values of critical temperature in the tran- 
sition ferromagnetic metals. Besides our results qualitatively reproduce the non-typical 
concentration dependence Curie temperature in systems Fei_ x Co a; S2 and Coi_a;Ni x S 2 [^ . 

In conclusion the correct taking into account of above mentioned matrix elements 
of electron-electron correlations allows to explain the behaviour of ferromagnetic transi- 
tion metals, their alloys and compounds both in the case of ground state and nonzero 
temperatures. 
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Figure captions 

Fig.l The critical values of zJ/w as a function of n at fixed value of U/w. Curves 1 
correspond to U/w = 1, T\ = t 2 = 0, curves 2 - to U/w — 1.2, r x = 0, r 2 = 0.15, curves 3 
. to U/w = 1.3, Ti = 0.15, r 2 = 0.1. 

Fig. 2 The ground state magnetization m GS as a function of n at [//«; = 1.2 and 
zJ/w = 0.1. Curve 1 corresponds to T\ = t 2 = 0, curve 2 - to Ti = 0, t 2 = 0.05 curve 3 - 
to T\ = 0.05, r 2 = 0, curve 4 - to Ti = r 2 = 0.05. 

Fig. 3 The Curie temperature as a function of n at values of U/w — 1 and r x = 0.05, 
r 2 = 0.1. Curves 1 corresponds to zJ/w = 0.4, curves 2 corresponds to zJ/w = 0.5, 
curves 3 corresponds to zJ/w = 0.6. 

Fig. 4 The dependencies of magnetization (solid curve) and Curie temperature (dashed 
curve) on electron concentration n at values of U/w = 1.5, T\ = 0, r 2 = 0.08 and 
zJ/w = 0.2 

Fig. 5 The Curie temperature as a function of U/w at half-filling. Curve 1 corresponds 
to zJ/w = 0.4, Ti = t 2 = 0; curve 2 corresponds to ^J/w = 0.4, ti = t 2 = 0.1; curve 3 
corresponds to zJ/w = 0.8, t\ = t 2 = 0; curve 4 corresponds to ^J/w = 0.8, r± = t 2 = 0.1. 

Fig. 6 The Curie temperature as a function of U/w at ^J/w = 0.4, T\ = 0.05, t 2 = 0.2 
and various band fillings. Curves 1, 2, 3, 4 and 5 correspond to n — 0.5, 0.75, 1, 1.25 and 
1.5 respectively. 

Fig. 7 The dependencies of magnetization and inverse magnetic susceptibility on re- 
duced temperature at U/w = 1, zJ/w = 0.5, T\ = 0.05 and t 2 = 0.1; the upper, middle 
and lower curves correspond to the case n — l,n — 0.5 and n = 0.2 respectively. 
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